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We propose a new method of quantization of a wide class of dynamical systems that originates 
directly from the equations of motion. The method is based on the correspondence between the 
classical and the quantum Poisson brackets, postulated by Dirac. This correspondence applied to 
open (non-Hamiltonian) systems allows one to point out the way of transition from the quantum 
description based on the Lindblad equation to the dynamical description of their classical analogs 
by the equations of motion and vice versa. As the examples of using of the method we describe 
the procedure of the quantization of three widely considered dynamical systems: 1) the harmonic 
oscillator with friction, 2) the oscillator with a nonlinear damping that simulates the process of the 
emergence of the limit cycle, and 3) the system of two periodic rotators with a weak interaction that 
synchronizes their oscillations. We discuss a possible application of the method for a description 
of quantum fluctuations in Josephson junctions with a strong damping and for the quantization of 
open magnetic systems with a dissipation and a pumping. 
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The main goal of this paper is to point out the close 
connection between the description of an open quantum 
system (OQS) in the Markovian approximation by the 
Lindblad equation and the description of its classical ana- 
log (CA), if it exists, by equations of motion presented 
in a special form. To settle this connection in the eas- 
iest way one can use the correspondence principle (CP) 
between the classical and the quantum Poisson brack- 
ets. The principle was postulated by Dirac in his famous 
book 0]. Having this connection in hands we propose a 
new method of quantization of non-Hamiltonian systems 
without memory. The procedure of using of the method 
consists in the following. Let us assume that the quantum 
system of interest has the CA with known equations of 
motion. We also imply that these equations can be pre- 
sented in a special form prescribed by the method - the 
form allowing the quantization (FAQ). In this case, on 
the base of the CP one can formulate an efficient recipe 
for finding the concrete form of the operators that give 
the quantum description of this system in the framework 
of the Lindblad equation. Below, we call such a pro- 
cedure the semiclassical quantization of the dynamical 
systems (without memory). We should emphasize that, 
to a certain extent, the method of quantization proposed 
is a heuristic one. Therefore, with reference to concrete 
examples, we demonstrate in detail the specifics of its 
application. As such examples we consider three popular 
dynamical systems: 1) the harmonic oscillator with a fric- 
tion; 2) the oscillator with a nonlinear damping that sim- 
ulates the mechanism of the emergence of the limit cycle 
near the bifurcation point; and 3) the system of two pe- 
riodic rotators with a weak interaction that synchronizes 
their oscillations. In our opinion, the latter example is 
in a direct relation with the description of the quantum 
noise in a Josephson junction with a strong damping. 
Besides, we give a simple generalization of the method 
proposed that allows one to apply it for the semiclassical 
quantization of open magnetic systems described by the 



equation of the Landau-Lifshitz type with a dissipation 
and under a pumping. Now, let us go to the details. 

It is well-known that quantum mechanics in its original 
form (the Schrodinger equation for the wave function ^' 
of the system: ihd'^/dt = H'if) provides the consistent 
description of the evolution of the state of a closed micro- 
scopic system. Such a description is extended automati- 
cally to the systems with a given dependence of their en- 
ergy on time: H — H[t). If a closed quantum system has 
the classical analog, the substitution of the semiclassical 
wave function ^' = Aevip^iSci/K) into the Schrodinger 
equation allows one to obtain the well-known Hamilton- 
Jacobi equation for the action function Sci{q,t) (in the 
limit h ^ Q), see, e.g. In that way one turns from the 
quantum to the classical description of a closed system. 
If the microscopic system S we are interested in is an 
open one, i.e. it interacts with an environment (a reser- 
voir R) that, as a rule, contains a large number of degrees 
of freedom, the quantum mechanics yields the following 
recipe for the description of the evolution of the state of 
S. One should consider the composite system SR that 
includes S and the reservoir R (that contains all exterior 
degrees of freedom with which the system S interacts). 
The evolution of the density matrix of the closed sys- 
tem PS nit) is described by the quantum Liouville - von 
Neumann equation: 

—^ = -j[Hsr,PsbI (1) 

where Hgji = Hg + Hn + Hint is the Hamiltonian of 
the closed system SR. In case when the solution of Eq. 

is known, the desired density matrix psit) is found 
by the reduction of psnit) over the variables that de- 
scribe the reservoir: psit) — trfl/5sj?Xi). As a rule, for 
real systems with interaction this general recipe is un- 
realizable. Therefore, most often, the procedure of re- 
duction is applied directly to the the equation JQ) and 
then one tries to obtain a closed equation for ps (t) using 
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approximations based on the presence of a small param- 
eter in a concrete physical system. On this way a large 
number of interesting and important results for quantum 
systems with dissipation were obtained (see and ref- 
erences therein). Nevertheless, such an approach is not 
a universal one, and in many cases it cannot be realized. 
Regarding the standard approach one should note that 
it looks rather " uneconomical" from the point of view of 
the experiment, since the description of an open system 
S requires the full information about the system SR. In 
connection with that, an alternative approach based on 
the quantum theory of open systems |J| becomes more 
and more popular during last years. In that approach, 
one attempts to describe the behavior of the OQS using 
only the information about the system S. The approach 
is based on the general restrictions imposed by the quan- 
tum mechanics on the evolution of psit): the linearity, 
the positivity, and the conservation of the trace of psit)- 
One of the most important results obtained in the frame- 
work of this approach is the Lindblad equation. This 
equation realizes the quantum scheme of the description 
of the evolution of the OQS in the Markovian approxi- 
mation. The Lindblad equation for the evolution of the 
density matrix of OQS psit) has the following form |^: 



where, by definition. 



N 



dps * rw - 1 

(2) 

The first term in the right hand side of Eq. Q de- 
scribes the internal (Hamilton) dynamics of S, while the 
other terms model its interaction with the reservoir R 
(TV is the number of channels of the interaction of S with 
R). We emphasize that the information determined in 
the standard approach by the structure of the reservoir 
and by the type of the interaction between S and R is 
included now in the operators Rj, R^ that acts in the 
space of states of the system S. 

As an example, let us consider the situation when N = 
1 and the operator R is the Hermitian one R = i?+ = O. 
It turns out that in this important special case Eq. ^ 
describes the decoherence of the quantum system inter- 
acting with the " meter" (macroscopic or mesoscopic one) 
that measures continuously the observable O, see 6j. 

Let us assume now that there exists a quantum sys- 
tem S (of the Markov type) which is described by Eq. 
^ with some /f, Rj and , and the system has the 
classical analog. The question emerges: how can one find 
the dynamics, i.e. the equations of motion, of that CA 
from the Lindblad equation for 5*? The answer on this 
question can be obtained with the help of the correspon- 
dence principle formulated by Dirac in According 
to Dirac, the quantum Poisson bracket (commutator) of 
two operators A and B is connected with the Poisson 
bracket {A, B} of their CA's - A(gr,Pr) and B{qr,Pr) by 
the relation 



lAm = V— — — (4) 

dqr dpr dpr dqr ' 



For the simplest operators A and B (for instance, A = qr 
and B = p^) the relation |(2Jl is exact one. In other cases 
it is fulfilled approximately, at least, in the first order in 
h. Let us use this CP for finding the classical analog of 
the Lindblad equation. We imply that in the limit ft — > 
the distribution function fs{<lr,Pr) in the phase space of 
the system S corresponds to the density matrix ps- We 
will consider the simplest case, when S has one degree of 
freedom and iV, the number of channels of the interaction 
of S and i?, is equal to one. Let the operators H , R and 
i?+ in Eq. Q be the functions of the operators q and 
p = {h/i)d/dq. Then, according to Eqs. l(2Jl and ||2Jl, 
one can write the following equation for the distribution 
function of CA of the system S: 



dfs 
dt 



{H,fs} + ih{Rfs,R} + ih{RJsR}, (5) 



where H{p,q), R{p,q) and R{p,q) are the classical 
analogs of H, R and J?"*", correspondingly. Using the def- 
inition and differentiating explicitly, we reduce Eq. 
(jSj to the following form: 

^ + Mq,p)^ + B{q,p)^ + c{q,p)fs = 0, (6) 



where the notations 



A{q,p) 
B{q,p) = 



dH 
dp 



ih\R 



OR 



dR 



dH 



ih\R 



dp 
dR 



^ dp 
j^dR 
dq 

C(q p) - 2ih ( _ dRdR 

^ dp dq dq dp 



(7) 



[A,B] = ih{A,B}, 



(3) 



are used. 

The first order partial differential equation © de- 
scribes the motion of the ensemble of systems (the phase 
liquid) in the space Ts = iq,p) of the system S. It is 
known (see 0) that the phase points (elements of the 
ensemble) move along the characteristics of Eq. (jS)) ac- 
cording to the equations of motion: 

p.Biq,p)^-- + ^n[^R--R-y (8) 

Now, let us compute the divergence of the phase ve- 
locity V = {q,p). By definition 

divv-^ + ^-2iftf— — — ^ (9) 
dq dp \ dq dp dq dp J 
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Using the relations ^ and ® we find that Eq. © is 
reduced to the standard form of the continuity equation 
of motion for the phase hquid in the space Ts 



^ + div(/sv) = 0. 



(10) 



Thus, with the help of the correspondence principle Q 
we arrive to the conclusion that in the limit h ^ the 
Lindblad equation Q turns to the continuity equation 
for the distribution function fs for the CA of the system 
S. The equations of motion for that CA (Eq. 0) follow 
directly from Eq. ||2Jl. This is the main result of the 
paper and we will use it for solving the inverse problem: 
the obtaining of quantum description of the dynamical 
system (without memory) from the equations of motion 
for its CA. 

It is convenient to present the equations of motion (jSJ 
in the form of one complex equation for the complex co- 
ordinate. Let us define the complex variables z and z*: 



iplg 
h 



iplg 
h 



V2 



V2 



(11) 



where Iq is an arbitrary parameter with the dimension of 
the length. It follows from the definition that for 
any choice of Iq the Poisson bracket of the dimensionless 
variables z and z* is 



According to the CP (see ®), [z, i+] = ih{z,z*} = 1. 
The latter relation allows one to identify the operators z 
and z+ with the standard Bose operators a and , that 
satisfy [a, a"*"] = 1. We note that in the case when the 
polynomial functions H , R and i? of 2: and z* contain the 
terms of the form z'^{z*y the ordering of the operators 
d and a+ should be fixed and necessary symmetrization 
should be fulfilled under the transition to the quantum 
analogs H, R and Below we demonstrate how this 
problem is resolved with reference to concrete examples. 
Note, that the same problem emerges under the quanti- 
zation of the Hamilton systems, as well. 

3. The operators H, R and i?+ should be substituted 
into the Lindblad equation ^ for the evolution of the 
density matrix. The correspondence principle guarantees 
that the equation obtained by this method will give a cor- 
rect quantum description of the evolution of S, at least, 
with the accuracy up to the first order in h. Therefore, 
the procedure described can be referred as the method 
of the semiclassical quantization of dynamical systems 
(without memory). We note that the most delicate point 
of the method is the presentation of the equations of mo- 
tion in the FAQ - it is difficult to formalize this point. 

The rest of the paper is devoted to the study of con- 
crete physical systems that can be quantized by the 
method proposed. 
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Using the system of equations of motion (jS]) for q and p, 
we find the desired dynamical equation for z: 



dz 



.dH 

dz* 



oz* oz* ) 



(13) 



It is implied here that the functions H{q,p), R(q,p) and 
R{q,p) in Eq. (|13|l are expressed as the functions on z 
and z* with the help of relations inverse to Eq. 



V2 



ih 



iz*-z). 



(14) 



The equation of motion for the complex coordinate z 
written in the form (|13|l is referred as the one presented 
in the form allowed the quantization (FAQ). Now, let 
us give the recipe (based on the CP) for the transition 
from the description of the CA of the system S by the 
equation of motion to the quantum description of S by 
the Lindblad equation. The recipe proposed consists of 
three consecutive steps. 

1. The input dynamical equations should be presented 
in the FAQ written for the complex coordinate z in the 
form H13|l . Such a representation determines automat- 
ically the functions H{z,z*), R{z,z*) and R{z,z*) en- 
tered in the FAQ. 

2. One finds the quantum analogs H, R and i?+ of 
that functions. For that one should replace the coordi- 
nates z and z* to the corresponding operators z and z~^ . 



Example A. Harmonic oscillator with friction 

This simplest example of a quantum dissipative sys- 
tem was considered more than once in the literature. In 
particular, in the recent review by Mensky ^ the equa- 
tion of the Lindblad type for the density matrix of such 
an oscillator was obtained by the method of restricted 
Feynman integrals (the method of quantum corridors). 
Let us show how simple the derivation of this equation 
(up to the notations) by the method proposed here. The 
equation of motion for the harmonic oscillator with fric- 
tion has the well-known form: 



mq + 2^q + muj'^q = 0, 



(15) 



where 7 is the friction coefficient, ujq, the oscillatory fre- 
quency, and m, its mass. Eq. (|15|l can be presented as 
the system of two first order equations: 



P 

9 = — , 
m 

27 2 

p — mojQq. 

m 



(16) 



Let us consider the complex coordinate z — {q/lQ + 
iplo/h)/\/2. In this problem it is convenient to choose 
the oscillatory length as the parameter Iq = h/mojo. 
At such Iq, using Eq. Ijltill . we obtain the equation of 
motion for z in the form: 



dz 
'dt 



-itooz — X{z — z*), 



(17) 
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where the notation A = 7/m is introduced. 

Let us show that Eq. (|17|l can be presented in the FAQ. 
We define two functions of the variables z and z*: the 
real one H = Hluo\z\'^ + ihX[{z*)'^ — z^]/2 and the complex 
one R = y/X{z cosh u — z* sinh u) , where u is an arbitrary 
parameter. Then we check directly that the right hand 
side of Eq. H18|l can be written in the form 

-..o.-A(.-z*) = ^- — -i?— j. (18) 

Therefore, in accordance with the method proposed, the 
quantum description of the harmonic oscillator with fric- 
tion is given by the Lindblad equation 

^ = -^[7?,p] + [7?p,i?+] + [i?,pi?+], (19) 

where the quantum analogs of the functions H{z, z*) and 
i?(z, z*) are determined as 

R = VX{acoshu — a^sinhw). (20) 

Equation IjlQI) together with the expressions for H and 
R (Eq. 1(201)) coincides with the result Eq. (51) of Ref. 
up to the notations. Using Eqs. ifTiHl and ((201), one 
can derive the equations of motion for the first {q, p) and 
the second (g^, p^, QP Pq) momenta of q and p - the 
basic physical quantities of the system S. But we will 
not do that here, since Ref. Q contains such a deriva- 
tion and the subsequent analysis (including a discussion 
on the physical meaning of the parameter u). We have 
considered the harmonic oscillator with only one aim: 
to illustrate all the specifics of the method of quantiza- 
tion on the example of the simplest model. Let us now 
turn to more interesting examples. To simplify the nota- 
tions, we will use the system of units with h = I. Also, 
the dimensionless coordinates x = q/lo and y = pIq/II 
will be used instead of the variables q and p. Conse- 
quently, the complex coordinate z in the FAQ is defined 
as z = (a; -|- iy)/V2. 

Example B. Quantum analog of the bifurcation of 
the emergence of the cycle 

Let us show how the quantum description of the oscilla- 
tor with a nonlinear damping can be given by the method 
proposed. The CA of such an oscillator simulates the 
mechanism of a smooth emergence of auto-oscillations 
near the bifurcation point. Similar bifurcation of the 
emergence of the cycle is typical for the situation when 
an equilibrium point losses its stability under a variation 
of the parameters of the dynamical system. In such a 
situation the growing up of small oscillations may switch 
the system to a new stationary state that corresponds to 
the closed trajectory (limit cycle). The following system 



of equations gives the mathematical model of the behav- 
ior of the dynamical system near the bifurcation of the 
emergence of the cycle (see 0): 

— = ujy + Xx - ^ix{x +y ), 
'^ = -ux + \y~^iy{x' + y^). (21) 

The system of the equations of motion Eq. H21f) can be 
written in the complex form: 

dz 

— ^-iuoz + Xz~2^z\z\^. (22) 

Equation H22|l can be presented in the FAQ. To show this 
we introduce the functions H = hujz*z, Ri — \/Xz* and 
i?2 — Then we check directly that the right hand 

side of Eq. H22I) can be rewritten in the form 

„ , ,2 dH 
— ibjz + Az — 2/iz|z| = 

According to the recipe of quantization described above, 
the Lindblad equation for the evolution of the density 
matrix of the quantum analog of the system (|21|l can be 
written in the form 

^ = + {Rxp,RX\ + {KpRX\ 

+ [^2/5,i?^] + [i?2,/5i?2+] (24) 

with ¥1 ~ bjo^a^ Ri ~ \/Aa+, and R2 — y^a^- 

From physical reasons it is clear that the motion of the 
system (I21f) along the limit cycle should correspond to the 
stationary solution of Eq. (|24|l pst, such that dpst/dt = 0. 
By analogy with the classical solution, for which jzp = 
const, one can expect that pst will commute with the 
particle number operator N = a^a. Therefore, we will 
seek the stationary solution of Eq. H24() in the form pst — 
Sn=o I I' where \n) are the eigcnfunctions of the 
operator N . Positive coefficients p„ satisfy the additional 
normalization condition X]ri=o Pn = ^- Substituting this 
expansion into Eq. (|24|l we obtain the system of recurrent 
relations for the coefficients pn'- 

2A[np„_i - {n + l)p„] + 
2/i[(n + 2)(n + l)p„+2 - n{n - l)p„] = 0. (25) 

It is convenient to introduce the generating function for 
the coefficients p„. By definition G{u) = J^'^^o Pn''^"' ■ 
One can check by the direct substitution that the system 
(|25|l is equivalent to the following second order differen- 
tial equation for the function G{u): 

{l+u)—-yu^-vG{u) = 0. (26) 
aw^ du 

In Eq. (|26|l the notation v — \/ p,\s used. 
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Let us transform Eq. H26|l to the canonical form. To 
do this we introduce the auxiliary variable v — ^{l + u). 
Then d/du — vd/dv, and Eq. 1)26(1 becomes 

Eq. (|27|l coincides with the equation for the confluent 
hypergeometric function y = $(a, c, x): 



x)^~av{x)^0 (28) 



(see ^3 ) if O'^s sets c—v and a = 1 . 

Therefore, the solution of Eq. 1)26(1 that satisfies all 
conditions of the problem can be written in the form 



G{u) 



$(l,t/,t/(l+M)) 



(29) 



In the solution Eq. ((29(1 we take into account that the 
normalization condition X]n=o = 1 is equivalent to the 
relation G(l) = 1 for the generating function G{u). 

Differentiating the expression Eq. ((29(1 one can find 
the average values of all physical quantities of interest in 
the stationary state. For instance, the average number 
of quanta n generated in the stationary state is given by 
the relation 



^ dG{u) 
2^ npn = 

n=0 



du 



(30) 



and the quantity n 



2 — /j-i - 



(n — l)n is connected with the 



second derivative of G{u) by the relation 



_ (PGiu) 



du^ 



(31) 



To characterize the properties of the distribution of the 
number of quanta |p„ | we will use the Mandel parameter 
Q defined as (see [lljj ) 



- 1 = 



n(n — 1) _ 
n. 



(32) 



The distributions for which Q > are called the super- 
Poisson ones, while the distribution with Q < 0, the 
sub-Poisson ones. The Poisson distribution corresponds 
to Q = 0. Using the properties of the confluent hyper- 
geometric function, one can easily find that at v = 1 the 
expression Eq. ((29(1 corresponds to the Poisson distribu- 
tion with n = 1, at 1/ > 1 it describes the super-Poisson 
distribution, and at < 1, the sub-Poisson distribution. 
We note also that for integer v the generating function 
G{u) can be presented through elementary functions. For 
instance, at = 1 it reads as G{u) = e"~^, aX v — 2, 
G{u) = (2/sinh2)(sinh(l w)/(l w))e"-i etc. Thus, 
the method proposed allows one to describe in detail the 
stationary state of OQS which classical analog simulates 
the bifurcation that occurs very frequently in the nonlin- 
ear mechanics. 



Example C. Quantum analog of the phase 
synchronization of two auto-oscillating rotators 

In this part of the paper we consider the specifics of ap- 
plication of the method for the problem of quantization 
of dynamical systems with several degrees of freedom. 
While such a generalization looks quite obvious, never- 
theless, additional interesting opportunities emerge, that 
will be analyzed in detail. As a distinctive example that 
is of interest by itself, we consider the dynamical system 
of two interacting auto-oscillating rotators. In the case 
when their frequencies (in the absence of the interaction) 
are rather close to each other the switching on the inter- 
action may result in the effect of synchronization, i.e. the 
locking of the frequencies and the phases in the bound 
system. The classical description of the effect of synchro- 
nization is well-known ■ Here we are interested in the 
possibility of its quantum description. We restrict the 
consideration by the case of weak interaction between 
the periodic oscillators, when the method of phase dy- 
namics (see can be used for their description. For 
the simplest form of the coupling between the oscillators 
this method yields the following system of equations of 
motion for their phases ipi and ip2 

dip 



—— =uji+a sm((^2 -Vi), 
dt 

— — = UJ2 + asm((^i - ip2), 
dt 



(33) 



where loi (0^2) is the frequency of motion the first (sec- 
ond) oscillator in its limit cycle, and a is the coupling 
constant that determines the strength of the interaction 
between the oscillators. 

It is implied in this approximation that each oscilla- 
tor (rotator) moves along the limit cycle with a constant 
amplitude, i.e., the additional restrictions for the system 
Eq. ((33(1 {dri/dt = dr2/dt — 0) are assumed. Taking 
into account these restrictions explicitly, one can write 
the system Eq. ((33(1 in the following equivalent form: 



uji + \r1r2 sin(((52 - 



dj£i 
dt 

= CJ2 + Arir2Sm((^i - (^2), 
dri _dr2 _^ 
dt ~ dt ^ ' 



(34) 



The system Eq. ((34(1 can be presented in the form of two 
dynamical equations for the complex coordinates zi — 
[xi + iyi)/V2 = ne^-^Vv^ and Z2 = {x2 + iy2)/^ 



r2e" 



VV2: 



dzi 
~dt 

dZ2 

~dt 



iUJlZl + \zi{zlz2 — Z2Z1), 
iuj2Z2 + \Z2{Z2ZI - Z\Z2)- 



(35) 



Indeed, if one substitutes the expressions for and Z2 
into Eq. 1(35(1 and then separates the real and the imagi- 
nary parts, one obtains just the system ((35() . Let us show 
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that Eq. (|35|l can be presented in the FAQ. To do this 
we choose the following functions H and R: 



H 



I |2 I |2 I ^'^l^lP / * * \ 

-UJl\zi\ H 1^ {z^Z2-Z2Zi) 



(Z2Z1 - Z^Z2), 



R 



^1 



^2 



Z^Zl - Z2Z^) 



(36) 



One can check directly that the right hand sides of the 
equations of the system are expressed as 

lUJlZl + \Zl{Z^Z2 - Z2Z1) = 

^dR _^dR_ ' 
dzl dzl 



iuJ2Z2 + \Z2{Z2ZI 



R 



ZIZ2) 

dR 

dz* ' 



dH 



R 



dR 

dz* 



(37) 



Then, following the method of quantization proposed, 
one can pass to the corresponding Lindblad equation in 
which the operators H and R are obtained from the 
functions 136|l by the substitution zi ^ ai, zl ^ a-i , 
Z2 ^ 0-2, Z2 ^ ■ But, there is another formulation of 
this problem, that is especially useful for the quantiza- 
tion of open magnetic systems. Let us describe in short 
this formulation. 

Let the dynamical systems with two degrees of freedom 
be determined by its equations of motion presented in 
FAQ with some functions H and R 



dzi 
It 
dz2 
~dt 



.dH 



dH 



R 



dR 



R 



dR 



dz 



R 



dR 

dl 



- R 



.dR 
d^ 



(38) 
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Instead of the complex coordinates zi and Z2 we intro- 
duce three real valued variables Ix = (^1^2 + Z2Zi)/2, 
ly = i{z^zi - zlz2)/2, and = {\zi\^ - \z2\'^)/2. We im- 
ply that the Poisson brackets for the coordinates zi, z 
Z2, and Z2 satisfy the relations {za, zp} = 0, {z*, z*p} — 
and {zq, z*p\ = —iSap, where a,l3 = 1, 2. Then it follows 
from the definition of Ix, ly and Iz that the Poisson brack- 
ets for the components of the vector 1 = {Ixjyjz) have 
the form {h^lk} = f-ijklk {i,j,k = 1,2,3). It allows one 
to identify the vector 1 with the angular momentum of 
the system considered. If H and R in Eq. H38|l can be 
presented as functions of Ix, ly and l^, then, substituting 
the variables and doing a simple algebra, we find that the 
system Eq. (I38II is equivalent to the following equation of 
motion for the vector 1: 



'dt 



SH 



iR 



, SR 



-I- c.c. 



(39) 



Eq. I|39|l can be considered as a generalization of the 
standard equation of motion for the magnetic momentum 



(Landau-Lifshitz equation) for the case of the momentum 
interacting with the environment. One can see from Eq. 
H39|l that such an equation is valid under the condition 
of conservation of the total magnetic momentum under 
the motion. Physically, it means that the interaction of 
the open system with the environment is invariant with 
respect to rotations of the system. The representation of 
the equations of motion for the magnetic momentum in 
the form H39|l with given H{i), R(l) and ^(1) allows one to 
fulfill the obvious transition to the quantum description 
on the base of the Lindblad equation. To do this one 
should find the operators H, R and R^ , the quantum 
analogs of the functions H, R and R, correspondingly. 
As before, the procedure of finding of that operators is 
determined by the correspondence principle: one should 
replace the variables Ix, ly and Iz in the expressions for 
H, R and R to the operators Ix, ly and Iz that satisfy 
the commutation relations [li,lj] = i^ijklk- Then, if nec- 
essary, the expressions obtained should be ordered and 
symmetrized. 

Let us now show how this variant of the quantization 
method works in the problem we are interested in - the 
problem on the synchronization of two interacting rota- 
tors. First of all, we note that the quantities H, R and 
R in Eq. H36fl are expressed as functions of Ix, ly and Iz 
as follows 



H 



-{UJI - U!2)lz - 2Xlylz, 



R = \^{lz 



til. 



(40) 



We omit the teiui — {lji+lu2){\zi\^ + \z2\'^) /2 in the expres- 
sion for H in Eq. H4U|I . because this term is the integral 
of motion of the equations H35|) . 

The quantum description of the system we are inter- 
ested in is given by the following Lindblad equation: 



[Rp, R+] + [R, pR^ 



(41) 



in which H and R are the quantum analogs of the func- 
tions H and R (Eq.lgOl)): 



H ^ -{uji - uj2)iz - A (jyiz + Liy 

R — VX (j^z " 'iiy 



(42) 



Using the Lindblad equation 141|) one can find the equa- 
tion of motion for the average value of an arbitrary ob- 
servable A: A = trpA. The direct calculation of dA/dt 
yields the following result: 



dA 
'dt 



-i[A,H]+ R+[A,R] + [R+,A]R. (43) 



Using Eq. I|43|l we obtain the system of equations of 
motion for the first order momenta Ix, ly and Iz'- 



dlx 
dt 



-2XL 



6L, 
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dly 

dt 



Xly 2,X(^lj;ly -f" lylx^ X 1 

-XL. 



(44) 



In Eq. (|44|l we introduce the notation S = coi — uj2 for 
the difference of the frequencies of two rotators. 

One can see from Eq. H44|l that the equations of motion 
for thejirst o rder mome nta contain the second order mo- 
menta l"^ and l^ly + lylx- Writing the equations of motion 
for the second order momenta 



dll 



dt_ 

dXl 
dt 



2X{ll-ll)-5{lylx+lJy), 



- 2X(lJy+ly Qly + 2 X^ ~ l^) , 



dt 



^ 8X11 ~3X{lyll+ Illy) 



-5X{lyl, + IJy) ~ 2XlJyl, + 25 [Py - l^) - Xly[45) 

we note that they contain the third order momenta, etc. 
Thus, further simphfications are needed for obtaining the 
solution of the quantum problem. Below, we consider 
only stationary states of the system (for which the left 
hand sides of Eqs. (|^ and are zero) and specify the 
case of complete synchronization <5 = 0. To obtain the 
closed system of equation for the first order momenta and 
the second order momenta we decouple the third order 
momenta applying the approximation frequently used in 
such problems (see 0): 



ABC K AB ■ C + A ■ BC + AC ■ B - 2A ■ B ■ C. 

As a final result, we obtain the following system of equa- 
tions for the first and the second order momenta 



Iz — Ixly — lylx 

'•y 2 ' ^ ^' 



4' 



Ixly ~^ Ix ' ly '^Ix ' ijy^ 



+ ill - ID = 0, 



3/y • iy 2 {iy 



0] i . J . 12 

^^x^y ''X I ''y ''x 



'2ly I Ix 



ly — 0' 



(46) 



It follows from the last equation of the system (|46|) that 
ly — 0. This relation is valid also for the CA. Let us 
remind that in the classical limit {h 0) ly ^ ly — 
rir2 sin((^2 ~ fi) — > 0, since for the case of complete syn- 
chronization the phase difference ip = ip2 — (pi 0. But, 
in difference with the CA, where l^ = r^r|sin^ (p is also 

approaches zero, in the quantum case lyj^O for the sta- 
tionary state, and this quantity can be the measure of the 



quantum noise of the system. Let us find the value of l^ 
using the system (I46II and the angular momentum conser- 
vation law il + Py + il = L2 = l{l + l) = {N/2){N/2 + 1), 

where N ~ a^di + 0^02 is the total number of the exci- 
tations under the bosonic description of the system. Let 

us notate x = ly. It follows from (|46|l that x = lx/2 and 
it satisfies the simple equation 

n 2 3a; 
The solution of Eq. JITJ) is 



N 



(47) 



/2 



1 



+ 2N + 



1/2 



N 



(48) 



In Eq. (|48ll we take into account that the semiclassical 
approximation requires the fulfilment of the inequality 
iV> 1. 

Let us summarize the main conclusions that concern 
the latter example. 

1. The quantum description of the non-Hamiltonian 
system consisting of two rotators which oscillations be 
synchronized by their mutual interaction is possible not 
only with the use of standard Bose operators, but also in 
the representation, where the components of the operator 
of the angular momentum 1 = (lx,ly,lz) be used as the 
generators. 

2. Under the condition of complete synchronization 
(lui = UJ2) the level of the quantum noise in this system 
is determined by the average value of the operator ly . 

3. In the stationary state ly ~ N/8, where N is the 
total number of the excitations under the bosonic de- 
scription of the system 

It is necessary to note that the modification of the 
quantization method used in the latter example is based 
on the Schwinger representation [T^ for the angular mo- 
mentum operator in terms of two Bose operators. 

To conclude this consideration, we discuss briefiy the 
possibility of the application of the results of the last 
section to the study of quantum fluctuations in Josephson 
junctions with a strong damping. Usually, the dynamics 
of such a junction is described by the resistive model on 
the base of the following well-known equation 



hip 
2^ 



Ic simp. 



(49) 



In Eq. (|49|l the following notations are used: p — 
fi — y^2 is the phase difference between the superconduct- 
ing banks of the junction, lex, the input current applied 
to the banks, R, the Ohmic resistance of the junction, Ic, 
its critical current. Comparing Eqs. H49|l and H'6'3\i we see 
that these equations coincide with each other up to the 
notations. The regime of the complete synchronization 
uji = UJ2 in the model of coupled rotators corresponds to 
the case of zero input current I^x in the Josephson junc- 
tion model. As is easily seen, the parameter of the quan- 
tum noise ly is proportional to the average square of the 
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bias voltage in the junction and this quantity can be mea- 
sured experimentaUy. It is also obvious that the number 
of bosonic excitations N — a^ai + 0^0,2 in the Joseph- 
son junction is equal to the number of Cooper pairs in its 
banks. In the framework of the Ginzburg-Landau theory 
this quantity is simply connected with the order param- 



eter of the superconductor [15| . The explicit form of the 
relation H48|l for the Josephson junction and a detailed 
analysis of the quantum fluctuations by the method pro- 
posed will be given in further publications. 

I would acknowledge L.A.Pastur for the discussion of 
the results of the paper and valuable comments. 
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